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ied in the framework of the nonlinear Schrodinger equation. It is found that 
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termination of catastrophic self-focusing of the pulse with anomalous disper- 
sion. It is also shown that the coupling between pulses can give rise to spatio- 
temporal splitting of the pulse propagating in anomalous dispersion regime, 
but it does not necessarily lead to catastrophic self-focusing of the pulse with 
normal dispersion. For the limiting case when the dispersive term of the pulse 
propagating in normal dispersion regime can be neglected an indication (based 
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was found earlier in media with saturation-type nonlinearity. 
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1 Introduction 

The propagation of a dispersive light pulse in a planar waveguide with positive, instantaneous 
Kerr-type nonlinearity can be described by the (2+l)-dimensional nonlinear Schrodinger equation 
(NSE) §: 
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where the parameters £, r, £ are as defined in appendix A. 

Equation (|l|) is valid only for pulses in the picosecond range; for shorter pulses additional 
terms, due to a higher-order dispersion, for example, should be included. The last term in equation 
([!]) describes Kerr-type nonlinearity; second and third terms are associated, respectively, with 
diffraction, which causes spreading of the pulse in space, and first-order group velocity dispersion, 
which leads to temporal broadening of the pulse. Parameter a, which can be either positive (for 
anomalous dispersion) or negative (for normal dispersion), is the dispersion-to-diffraction ratio 
||. The spatio-temporal dynamics of the pulse depends, to a high degree, on the sign of this 
parameter. 

It is known that some solutions of the (2+l)-dimensional NSE can develop into a singularity 
of the electric field in the self-focus point. This phenomenon, known as catastrophic self-focusing, 
occurs simultaneously in space and time for pulses propagating in planar waveguides with ano- 
malous group velocity dispersion (equation (|l]) with a > 0) |ij ||, and also for dispersionless 
beams propagating in self-focusing bulk media (equation ([l]) with the dispersive term replaced 
by a diffraction term) Q when parameters of the system are above the threshold of catastrophic 
self- focusing Q , which is usulally computed with the aid of the method of moments ||, [?], |§| , the 
variational method || [)), and also numerical simulations Jl|, [lO], The occuurence of catas- 
trophic self-focusing is not only non-physical, it also prevents examination of the pulse behavior 
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behind the self-focus, for it emerges just as an artifact of approximations made when deriving 
the NSE. In order to avoid this limitation, either some nonlinear stabilization mechanisms such 
as saturation fl3| or non- locality jl3| of nonlinearity, Raman scattering [p3, plasma formation 
|l5| ], multiphoton ionization fllBfl , higher-order group velocity dispersion terms jlTj], an adequate 
composition of the above mentioned effects |18|, or a non-paraxial treatment of the process of 
self-focusing Jlj], |J should be included into consideration. However the standard 

paraxial NSE can still serve as the model equation for sclf-focusing in the case when parameters 
of the system are below the threshold of catastrophic sclf-focusing or, in the reverse case, for 
studying dynamics of a pulse/beam in the prefocal region. 

Another situation occurs when the pulse propagates in a normal dispersion regime. In this 
case the terms describing dispersion and diffraction have different signs and two different effects, 
spatial self-focusing and temporal self-defocusing, simultaneously influencing the propagation 
of the pulse. This causes the situation where, in the solution of the NSE (equation (1) with 
a < 0) neither singularity |Q nor localized steady-states occurs p5| . Moreover, this solution is 
accompanied by a breaking of spatio-temporal symmetry and a uniform structure of the pulse 
and can finally lead to an occurrence of several humps in the field distribution |25f| , splitting 
of the pulses into two sub-pulses |^6) , or splitting into several sub-pulses [^7| . It has also been 
reported that the presence of even very small normal dispersion can lead to the destruction of 
soliton breathers propagating in nonlinear planar waveguides pq . In the case of the (3+1)- 
dimensional NSE splitting of a pulse into two sub-pulses has also been observed 
while splitting into several sub-pulses predicted theoretically in p0| , |35| , [34], |35|, |36j has been 
confirmed experimantally by the authors of EG, p?L 
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Thus, depending on the sign of dispersion, a dispersive pulse propagating in a Kerr-type 
planar waveguide reveals different behaviour. Catastrophic self-focusing (in the framework of the 
NSE) takes place in the case of anomalous dispersion. For normal dispersion the typical process 
is spatio-temporal splitting. It seems interesting to study an interaction between two pulses co- 
propagating in such a medium, i.e. a Kerr-type planar waveguide, under the assumption that one 
of them propagates in a normal dispersion regime and another is in an anomalous regime. To 
the author's knowledge this problem has not been studied in the literature and the main purpose 
of this paper is to consider it. Note that interaction of spatially separated light beams whose 
evolution is modeled by a set of n (n > 2) nonlincarly coupled NSEs was studied by several 
authors |5j| [Io[ fyj 42 . Moreover, the importance of the interaction between two pulses in a 



nonlinear medium has been pointed out already by Agrawal in ]43| , where an intriguing effect of 
an induced focusing of two beams co-propagating in a self-defocusing medium has been reported. 

It is also known that neither for anomalous dispersion f44|| nor for normal dispersion |p5| 
do stable soliton-like solutions of the (2+l)-dimensional (and also (3+ l)-dimensional) NSE exist. 
This statement also concerns experimental results, since no soliton-like solution has been observed 
in pure Kerr-like nonlinear media with two or three transverse dimensions. From the point of 
view of applications, i.g. as elements of optical switching devices |4j|, the existence of stable 
soliton-like solutions is very important. Therefore, solutions to this problem has been already 
proposed by several authors: for example, it has been shown that soliton-like structures can 
be realized in media with saturation-type nonlinearity [12, 46, WA |48|], in photorefractive media 
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[|49|, gGj, in media with quadratic nonlinearity 51, |52 
^(2) _ ^(3) nonlinearity El], pq], and also in the limiting case of the discrete-continuous NSE 



which can model propagation of short optical pulses in an array of linearly coupled optical fibers 
p9| . In this paper we will consider another possibility of obtaining a self-trapped solution in two 
transverse dimensions, namely in a configuration of the (l-l-l)-dimensional NSE coupled to the 
(2+l)-dimensional NSE. 

We proceed as follows. In section 2, two coupled NSEs describing the co-propagation of 
two dispersive pulses in a nonlinear planar waveguide and basic equations following from the 
variational method will be introduced. Next, in section 3, the problem of catastrophic self- 
focusing will be considered. First, the influence of the parameters of the pulse propagating in a 
normal dispersion regime on the threshold of catastrophic self-focusing of the pulse propagating 
in an anomalous dispersion regime will be studied. We will also examine whether catastrophic 
sclf-focusing of the pulse propagating in a normal dispersion regime can occur as a result of the 
nonlinear coupling between two pulses. In section 4, which is devoted to the problem of spatio- 
temporal splitting, we will investigate whether the influence of the pulse propagating in a normal 
dispersion regime can enforce spatio-temporal splitting of the pulse with anomalous dispersion. 
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In the last section, section 5, we will focus on the limiting case when the dispersive term of the 
normal pulse can be neglected. In this case the problem of two coupled (2+l)-dimensional NSEs 
will be reduced to the system of a (l+l)-dimensional NSE coupled to a (2+l)-dimensional NSE. 
The main reason to study this configuration is to investigate a possibility of a stable, self-trapped 
solution. 

The interaction between pulses will be assumed to be limited to cross-phase modulation, a 
nonlinear effect through which the phase of an optical beam/pulse is affected by another propagat- 
ing beam/pulse and which can cause a redistribution of energy within each beam/pulse. Another 
effect, four-wave mixing, will be neglected, so that no energy transfer between both pulses will 
be taken into consideration. The analysis presented in this paper is based on the variational 
method [|60| and numerical simulations using the split-step spectral method J(HJ . Throughout the 
paper the pulse propagating in an anomalous (normal) dispersion regime will be referred to as 
the anomalous (normal) pulse. 



2 Basic equations 

The co-propagation of two optical pulses in a nonlinear planar waveguide can be described by 
two coupled nonlinear Schrodinger equations: 
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where the last terms represent cross-phase modulation, a nonlinear effect which causes a coupling 
between pulses and the terms before the last ones describe self-phase modulation. 

It is assumed that the subscript j = 1 (j = 2) denotes the anomalous (normal) pulse, hence 
o\ > and 02 < 0. The notations in equations (§a,b) are explained in appendix A. The initial 
conditions will be taken in the form of the Gaussian pulses 



*i(C = 0,r,£) = JK~exp 
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where C T j (Cgj) is the temporal (spatial) chirp of the j-th pulse, j = 1, 2. The parameter Kj will 
be called here the strength of nonlinearity of the j-th pulse (see explanation in appendix A) . 



2.1 Variational method 



It is known that the set of NSEs (equation 
given by 



,b)) can be obtained from the Lagrangian density 
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Following the variational method ||60|| let us choose a proper multi-parametric trial function for 
the solution of equation (||a,b). Since in this paper we consider the Gaussian initial condition 
(equation (jsh) it is natural to take as the trial function the Gaussian function: 
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(5) 



with 12 parameters: the complex conjugate amplitudes, Aj,A*, the temporal and the spatial 
widths, w T j,W£j, and the temporal and the spatial chirps, C T j,C^j, where j = 1,2. From the 
initial condition (equation (^)) it follows that Aj(( = 0) = ^/KJ, w T j(C = 0) = w^j(C = 0) = 1. 
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The evolution equations for the parameters of the trial function are obtained by varying the 
reduced Lagrangian 



(L) := J Ld£_dr, 



into which the trial function (equation (^|)) is inserted, with respect to the parameters of the trial 
function, Aj,A*, w T j, W^j, C T j, C^j. We obtain the following 12 coupled ordinary differential 
equations: 
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where X, := w T j(C)w^(C)|A(C)| 2 = Kj = const. From equations (6a) and (6b), which are actually 
the energy conservation laws for two pulses, Nj ■= J J \^j\ 2 drd^ — irTj, it follows that there is 
no energy transfer between the pulses. 

The set of equations (6a)- (9b) is rather complicated and only in the special case when g\ = 
a 2 = 1 , X 2 = is the analytical solution 
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available J62] . More general situations should be treated numerically, e.g. using the Runge-Kutta 
method |63[ |. Still, equations (7a)-(7d) can be simplified to one evolution equation. To proceed, 
let us first rewrite equations (7a)- (7d) in the form 
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where the potentials Vi(w T i,W£i,w T 2,'W£2) and V2(w T i,w^i,w T 2,w^ 2 ) read as 
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is a constant of motion. 

Again, using equations (7a)-(7d) it can be calculated that 
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where w := rX\ 



(here we assume a% 7^ 0, a 2 ^ 0, fx ^ 0). 



From equation (|l0[) one can easily get the evolution equation for w 
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3 Catastrophic self-focusing 

This section is devoted to the problem of catastrophic self-focusing, which can occur in the solution 
of the set of equation (||a,b). Our analysis is based on the variational method and numerical 
simulations and the comparison of the results of both. Note that once we have specified the 
threshold of catastrophic self-focusing we then know for which parameters of the system the NSE 
is valid and we can use this information in further research. 

From the point of view of analytical estimations, which can be done using the method of 
moments [|[ Q or the variational method |}5| catastrophic self-focusing is identified with a 
development of a singularity in the solution at a finite distance of propagation. 

Let us briefly discuss the case of a single pulse, i.e., let us make the assumption that I2 = 
and X := X\, a :— a\,w T :— w T i,W£ := w^i. Then we get that 
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whereas equations (|10|) and (y_l|) still remain valid. From equation ( |ll| ) it follows that the quantity 
w goes to zero on a finite distance of propagation when one of the following conditions is satisfied 



r w < 0, 
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(12) 



< -2y/Ww(0). 



A vanishing of w can be associated with a singularity of the solution of the NSE (equation 
([!])) only when dispersion is anomalous, a > 0, since only in this case the quantity w can 



5 



be interpreted as an average width of the pulse and the condition w — is equivalent to a 
simultaneous vanishing of both widths of the pulse. Therefore, for the Gaussian initial condition 
(equation (||) with k := «i) without initial chirp, C T (0) := C T i(0) = 0, Q(0) := Qi(0) = 0, i.e. 

= 0, catastrophic sclf-focusing of the pulse with anomalous dispersion will arise when 



for f£ 

the condition W < is satisfied, i.e. when 



K > KVcat = O + 1. (13) 

Note that the condition given by equation dl3| ) agrees with the results obtained in Q with the 
aid of the method of moments for an elliptic Gaussian beam. 

Another situation occurs in the case of normal dispersion, namely a vanishing of the quantity w 
means only that = erw|, therefore, nothing about catastrophic self-focusing can be concluded 
from equation (pi]). However, based on the method of moments it has been demonstrated that 
catastrophic self-focusing in this case does not occur p4| . 

In our numerical simulations catastrophic self-focusing is identified with a discontinuity of 
the phase </>(t, £, £) of the amplitude \1/ := l^le^ in the central point of the coordinate system, 
r = 0,£ = 0, and with non-monotonic behavior of the intensity l^l 2 in the central point after 
catastrophic self-focusing has been reached . The threshold of catastrophic self-focusing given 
by the numerical analysis |h], [m], |ll[ 

KNcat ~ <* + 0.885 

is lower than the one given by analytical estimations. 

Let us examine now two coupled NSEs given by equations (§a), (§b). We can consider three 
different cases: (i) both pulses propagate in an anomalous dispersion regime, o\ > 0, a 2 > 0; 
(ii) both pulses propagate in normal dispersion regime, a\ < 0, a 2 < 0; (iii) pulses propagate in 
different dispersion regimes, anomalous and normal, a± > 0, a 2 < 0. 

In the first case, when both pulses propagate in anomalous dispersion regimes the threshold 
of catastrophic self-focusing can be calculated in a similar way as it was done for a single pulse 
and is given by equation (12). For the Gaussian initial condition (equation (|j|)) without initial 

chirp, C r i(0) = Cfi(0) = C T 2(0) — C^ 2 (0) — 0, i.e. when ^ =0, catastrophic self-focusing 

occurs when the condition W < 0, which reads as 

rXi (01 + 1 - Ki - 2k 2 ) + l 2 {a 2 + ^i- rn 2 - 2mi) < 0, (14) 

is satisfied. 

Since vanishing of the quantity w, which can be interpreted as an average width of the pulses, 
is associated with a simultaneous vanishing of both widths of both pulses, then it can be concluded 
that when catastrophic self-focusing of one of the pulse occurs, it also occurs for the second one. 
This conclusion and the condition (pT ) written for the symmetric case, o\ = a-i = n = r = 1, agree 
with the results obtained in |59[ |4C|, 4l|, [l2| for two cylindrically symmetric, spatially separated 
beams whose distance vanishes. 

In the case of two pulses propagating in a normal dispersion regime the situation is simple: 
catastrophic sclf-focusing does not develop, even for very large strengths of nonlinearity of the 
pulses. 

The situation in more complicated when the pulses propagate in different dispersion regimes, 
anomalous and normal: the threshold of catastrophic self-focusing cannot be calculated from 
equation (pi]), therefore numerical solutions of equations (7a)-(7d) should be performed in order 
to analyse this problem. The first goal of our study is to examine an influence of the parameters 
of the normal pulse on the threshold of catastrophic self-focusing of the anomalous pulse. The 
parameters of the anomalous pulse have, therefore, been chosen in such a way that the relations 
Ki > Kvcat = l + ci (in the variational method), and k\ > KNcat ~ 0.885 + a\ (in the numerical 
simulations) are satisfied, which mean that catastrophic self-focusing of the anomalous pulse will 
develop when there is no coupling between pulses. Then the parameters of the normal pulse, i.e. 
the strength of nonlinearity, k 2 , and the dispersion-to-diffraction ratio, a 2 , are varied. We found 
that catastrophic self-focusing of the pulse propagating in an anomalous dispersion regime can 
be arrested by the influence of the pulse propagating in a normal dispersion regime. 
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The results following from the variational method are shown in figure |. The shaded area 
denotes the range of the parameters of the normal pulse, K2 and 02, for which catastrophic self- 
focusing of the anomalous pulse does not occur. It is evident that for small nonlinearity of the 
normal pulse, K2, the term describing cross-phase modulation of the anomalous pulse is negligible 
as compared with self-phase modulation. Therefore, the process of catastrophic self-focusing 
cannot be stopped and it takes place for all values of 02. When the strength of nonlinearity 
K2 increases, the influence of the normal pulse on the anomalous pulse through the cross-phase 
modulation term increases and then it is possible, for some values of the dispersion-to-diffraction 
ratio, \aviow(K2)\ < | cr 2 1 < |0Vu PP (k2)|, to stop catastrophic self-focusing. The lower threshold, 
\cviow( K 2)\, m the beginning decreases with an increase of the strength of nonlinearity of the 
normal pulse, «2- For sufficiently large nonlinearity, k 2 > Kyiow, the lower threshold becomes 
zero. The upper threshold, |oy U pp(/J2)|) increases with an increase of nonlinearity. The existence 
of the lower threshold can be explained as follows: when \a%\ < |ov/ OUI | is small, the dispersive 
term of the normal pulse is negligible as compared with diffraction. Therefore, the most important 
role in the propagation of the normal pulse is played by self-focusing ,which not only does not 
lead to an arresting of catastrophic self-focusing, but even additionally enhances it. A similar 
situation is known, for example, in a configuration of two beams, which co-propagate in a bulk 
medium and have the same amplitudes fl39| : namely the critical value of nonlinearity necessary for 
catastrophic self-focusing is three times smaller than in the case when they propagate as single 
pulses, as can be calculated i.g. from equation (|l4|). On the other hand, while for large | cr 2 1 , 
we have a broadening of the normal pulse with a significant spreading of energy out from the 
centre of the coordinate system, £ — 0,t = 0, while for the anomalous pulse there is a tendency 
of energy to concentrate in the centre. Then the overlap of two pulses becomes negligible, so 
that the coupling between them through cross-phase modulation is very small and catastrophic 
self-focusing of the anomalous pulse can not be stopped by the influence of the normal pulse. 

The results obtained with the aid of the numerical calculations are shown in figure ^. They 
confirm predictions of the variational method. Namely, catastrophic self-focusing of the anoma- 
lous pulse can be arrested by the pulses propagating in a normal dispersion regime when the 
strength of nonlinearity is sufficiently large, K2 > kniow, and the dispersion-to-diffraction ratio 
satisfies the relation \a N i ow (K 2 )\ < |tr 2 1 < \<Jnu PP {k2)\- 

Another question is whether the nonlinear coupling between pulses can cause catastrophic self- 
focusing of the pulse propagating in a normal dispersion regime. As it has been already mentioned, 
in the case of two simultaneously propagating pulses with anomalous dispersion catastrophic self- 
focusing of one of the pulse is associated with catastrophic self-focusing of the other one, and 
both widths of both pulses go to zero simultaneously when catastrophic self-focusing occurs. 

However, the variational method demonstrates that in the case discussed here catastrophic 
self-focusing of the anomalous pulse does not necessarily lead to catastrophic self-focusing of the 
normal pulse. Namely, when catastrophic self-focusing of the anomalous pulse occurs, the normal 
pulse can demonstrate, depending on the parameters of the system, two different characteris- 
tics: (i) both widths of the pulse initially decrease reaching a minimum on a certain distance of 
propagation and then they start to increase; (ii) the spatial width of the pulse vanishes to zero 
on a finite distance of propagation whereas the temporal width initially decreases, reaching a 
minimum on a certain distance of propagation, and then it increases. In particular, the case (i) 
can be realized for the following parameters of the system K\ = 3, K2 = 3, a\ = 1, 02 = —7, while 
the case (ii) occurs, for example, when k\ = 3, K2 = 3,cti = 1,02 = —1. An effect, similar to 
(ii) has been also observed in the case of a pulse that propagates in a bulk medium with normal 
dispersion and whose dynamics is modeled by the (3+l)-dimensional NSE: namely catastrophic 
self-focusing of this pulse occurs when only the spatial widths vanish to zero while the temporal 
width is not allowed to reach zero value [[56] . 

The numerical simulations have not confirmed the results of the variational method concerning 
a possibility of catastrophic self- focusing of the normal pulse whose spatial width vanishes to zero 
on a cartain distance of propagation, £, and whose temporal width is left larger than zero. 
However, no definite statement that this effect is prohibited can be made either. Additional 
calculations should be performed to clarify this question. 

We can therefore conclude, based on the variational method and the numerical simulations, 
that catastrophic self-focusing of the anomalous pulse does not necessarily lead to catastrophic 
self-focusing of the normal pulse. 
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Figure 1: The results of the variational method displaying the dependence of the threshold of 
catastrophic self-focusing of the pulse propagating in an anomalous dispersion regime, ^i, on the 
parameters of the pulse propagating in a normal dispersion regime, ^2- The shaded area denotes 
the range of the parameters, the strength of nonlinearity, K2, and the dispersion-to-diffraction 
ratio, cr 2 , for which catastrophic self- focusing does not occur. The parameters of the anomalous 
pulse have been chosen in such a way that they are above the threshold of catastrophic self- focusing 
in a single propagation regime: i.e. for (a) n x = 2.2, a x = 1.0, for (b) Ki = 2.3, a x = 1.0. 
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Figure 2: The results of the numerical simulations displaying the dependence of the threshold 
of catastrophic self- focusing of the pulse propagating in an anomalous dispersion regime, ^i, on 
the parameters of the pulse propagating in a normal dispersion regime, ^2, i-e. the strength 
of nonlinearity, K2, and the dispersion-to-diffraction ratio &2- Full circle (empty triangle) points 
denote the occurrence (lack) of catastrophic self-focusing. The parameters of the anomalous pulse 
have been chosen in such a way that they are above the threshold of catastrophic self-focusing in 
a single propagation regime, i.e. «i = 2.0, o\ = 1.0. 



4 Spatio-temporal splitting 

In this section the problem of spatio-temporal splitting is discussed in more detail. The origin of 
spatio-temporal splitting of a single pulse propagating in a normal dispersion regime in Kerr-type 
planar waveguides p7| |G7], [|(J or bulk media J2{| |3(], [u], [32| is the fact that spatial self-focusing (in 
one or two dimensions) and temporal self-defocusing act simultaneously during the propagation. 
Therefore, in space, there is a tendency of energy to concentrate in the centre of the coordinate 
system, r = 0, £ = 0, whereas in time a spreading of energy away from the centre takes place. 
When both effects are combined, local focusing areas develop away from the centre and, as a result, 
spatio-temporal splitting of the pulse into several sub-pulses takes place p7f . The number of sub- 
pulses emerging in this way, for a sufficiently large propagation distance, is proportional, as it has 
been proposed in [^7|, to the order of the temporal soliton, N := y/ n/a, where the parameters 
of the (2+l)-dimensional NSE (equation ([[])) a and n denote respectivelly the dispersion-to- 
diffraction ratio and the strength of nonlinearity. Specifically, splitting of the pulse into two 
sub-pulses has been observed |26| in the system with parameters C = 2, k = 4, a = —0.1, while in 
p7[ splitting of the pulse into three sub-pulses has been obtained for ( = 0.15, k = 100, a = —3. 

Although we have not verified the statement that the number of sub-pulses is proportional to 
the order of a temporal soliton, N, since it was not the purpose of our study, we have observed 
that in the case of a single pulse with normal dispersion the tendency of the pulse to split increases 
when the strength of nonlinearity, k, increases and when the dispersion-to-diffraction ratio, a, 
decreases. Still, it remains for us an open question whether the humps in the filed distribution of 
a pulse mentioned in p5[ can be identified with sub-pulses whose existence is demonstrated in this 



paper and in 27, 26]] . If this is not the case, we could take the opportunity to speculate that the 
splitting of a pulse could not be observed by the authors of |25| since the strength of nonlinearity 
used by them was relatively weak, k ~ 1.76 (while the dispersion-to-diffraction ration was chosen 
to be a w —0.32) and only small local humps, instead of full pulse splitting, could be detected. 

Some agreement between numerical and variational solutions of the NSE with normal dis- 
persion has been demonstrated in the literature. For example, in [^H p7| it has been shown, 
respectively, that in both methods the number of oscillations of the peak amplitude of the pulse 
is the same and there is a similarity in evolution of the average square widths of the pulses. 
Nevertheless, in all cases when spatio-temporal splitting of pulses has been observed, numerical 
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simulations have been used p7L pOl [32l [37l 67 1. The variational method is not appropriate to 



predict the splitting of pulses, since it requires the solution to have a shape which does not 
change in propagation. When a Gaussian function is chosen as the initial condition, as we have 
done in this paper, it is difficult (if not impossible) to guess a trial function function which would 
satisfy the initial condition and also could describe spatio-temporal splitting of the pulse. Here, 
it is worthwhile to recall that the variational method cannot be applied to predict, for example, 
the formation of higher-order solitons in planar waveguides or optical fibers [poj . 

Since the variational method is not applicable to the study of spatio-temporal splitting of two 
pulses propagating simultaneously in a nonlinear planar waveguide, the results of this section 
are due to numerical simulations. Figures |3| and || show the spatio-temporal dependences of 
intensities of both pulses, the anomalous one (figure ||(a) and [|(a)), and the normal one (figure 
||(b) and ||(b)), for different longitudinal variables, (. Parameters of the pulses were chosen in 
such a way that when they propagate as single pulses the following effects take place on large 
propagation distances: (i) symmetric, spatio-temporal broadening of the anomalous pulse without 
occurence of catastrophic self- focusing (see figures ^(c) and [|(c)); and (ii) large asymmetrical, 
spatio-temporal broadening of the normal pulse without splitting into sub-pulses (see figures 
||(d) and [|(d)). The case (i) occurs when the conditions o\ = 1 and k,\ < Kycat = 1 + Ci 
are satisfied, while the case (ii) docs when the strength of nonlinearity, K2, is sufficiently small. 
When the pulses propagate simultaneously, i.e. there is a nonlinear coupling between them, the 
situation becomes qualitatively different, as it can be seen from figures |(a) and |(b). Namely, 
spatio-temporal splitting of both pulses can develop, so that for the propagation distance £ = 2 
the anomalous (normal) pulse becomes divided into n ~ 3 (n > 10) sub-pulses. The effect of 
splitting of the anomalous pulse, which does not occur when it propagates as a single pulse, 
can be explained as follows. When the nonlinear coupling between pulses through cross-phase 
modulation is present one pulse can induce a redistribution of energy of the other pulse. Therefore, 
if there are some local focusing areas in the distribution of energy of one pulse, energy of the other 
pulse tends to concentrate there. Such a tendency has been already pointed out by Agrawal, who 
has observed the occurrence of local focusing areas in the distribution of energy of two beams 
which co-propagate in a defocusing nonlinear medium 0] . 



5 The limiting case of vanishing dispersion of the normal 
pulse 

In this section we consider the limiting case when the dispersive term of the normal pulse can be 
neglected. We apply the variational method and numerical simulations and compare their results. 
We assume that the initial condition has the shape of the Gaussian function given by equation ([|) 
and concentrate basically on the question as to whether there exists a stable self-trapped solution 
of the above-mentioned system of equations. 

First we briefly discuss the case when the pulses propagate in a planar waveguide separately, 
i.e. when there is no coupling between them. Specifically, we consider (i) the propagation of a 
pulse with anomalous dispersion and (ii) the propagation of a dispersion-less beam. The case 
(i) can be described by the (2+l)-dimensional NSE, which does not have stable, self-trapped 
solutions. Thus, depending on parameters of the system, either spatio-temporal spreading of the 
pulse or catastrophic self-focusing develops. The case (ii) is described by the (l+l)-dimensional 
NSE which depends only on one transverse variable £ and being an integrable system possesses the 
familiar soliton solution given by sech function Q. Taking the Gaussian function (equation (^|)) 
which depends on two transverse variables r and £ from the variational method we obtain that 
the temporal width of the pulse is constant while the spatial width oscillates. These oscillations 
are due to the fact that the shape of the Gaussian trial function differs from the exact soliton 
solution given by the sech function [6q] . However, numerical simulations lead to a slightly different 
behavior. Namely, the temporal width of the pulse appears to oscillate in sinchronization with 
the spatial width. Amplitudes of both oscillations decrease with the longitudinal variable C and 
vanish at finite C when the spatial soliton is formed [j69| . 

Now, let us take into account the nonlinear coupling between pulses and discuss briefly the 
aspect of catastrophic self-focusing. Based on the variational method we have observed that when 
parameters of the anomalous pulse are chosen in such a way that catastrophic self-focusing does 
not occur when it propagates as a single pulse, i.e. when the condition k\ < Kvcat = 1 + &i is 
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C = i-o 

Figure 3: The spatio-temporal dependence of the intensity of the anomalous pulse (a) 
co-propagating with the normal pulse (b) and the spatio-temporal dependence of the intensi- 
ties of both pulses, anomalous one (c) and normal one (d) when they propagate separately. The 
following parameters of the system have been chosen: o\ = 1.0, Ki = 1.88, oi = — 0.1, k 2 = 2.0; 
the propagation distance C = 1-0. 
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Figure 4: The spatio-temporal dependence of the intensity of the anomalous pulse (a) 
co-propagating with the normal pulse (b) and the spatio-temporal dependence of the inten- 
sities of both pulses, anomalous one (c) and normal one (d) when they propagate sepa- 
rately. The parameters of the system are the same as in the case considered in figure 3: i.e., 
<7i = 1.0, K\ = 1.88, <T2 = —0.1, k 2 = 2.0; the only difference is that the distance of propagation is 
larger, ( = 2.0. 
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satisfied, then also in the case of the anomalous pulse coupled to the normal pulse no catastrophic 
self-focusing of both of them occurs, even for very large strength of nonlinearity of the normal 
pulse, i.g., k-2 = 20. However, when the condition k\ < Kvcat = 1 + ci is not fulfilled a similarity 
to the case discussed in section 3 can be found, namely three different behaviours of the pulses 
can be observed: (i) no catastrophic self-focusing of both pulses; (ii) catastrophic self-focusing 
of the anomalous pulse; (iii) catastrophic self-focusing of both pulses, the spatial width of the 
normal pulse vanishes to zero while his temporal width remains larger than zero. 

The above results obtained in the variational method have been verified in the numerical 
simulations except for the case (iii): that is a development of catastrophic self-focusing of the 
normal pulse whose spatial width venishes to zero on a certain distance of propagation, £, and 
whose temporal width is left larger than zero has not been observed, but also no definite statement 
that this effect is prohibited can be made. Rather, we are interested in possibilities of formation 
of self-trapped solutions and, therefore, we will restrict our analysis to parameters of the pulses 
which assure that no catastrophic self-focusing occurs, i.e., that condition k% < KNcat ~ 0.885+<Ti 
is satisfied. 

From the variational method it follows that the evolution of the normal pulse coupled to the 
anomalous one is essentially similar to that of the single normal pulse. Namely, the temporal 
width of the pulse does not depend on the longitudinal variable, £> as is seen from equation (8c) 
with the neglected dispersion of the normal pulse 02 — 0, while the spatial width of the pulse 
undergoes periodic oscillations (see figure |^(b)). The propagation of the anomalous pulse coupled 
to the normal one is, however, qualitatively different to the behaviour of a single anomalous 
pulse. Namely, both temporal and spatial widths of the pulse undergo periodic oscillations (see 
figure |^(a)). Therefore, neither spatio-temporal spreading nor catastrophic self-focusing of the 
anomalous pulse can develop and a self-trapped solution arises. Note that a similar self-trapped 
solution was found in the case of (2+l)-dimensional NSE with the saturation of nonlinearity p2|. 

We also performed numerical simulations for the case of simultaneously propagating pulses. 
The results are displayed in figure ^| from which it is evident that the temporal and spatial widths 
of both pulses oscillate in sinchronization, with the amplitude of the temporal oscillations smaller 
than the amplitude of the spatial ones. Unfortunately, the numerical calculations are rather 
labourious and we have not yet been able to calculate evolution for longer longitudinal variables, 
£ > 2, so that we do not know whether the amplitude of oscillations decreases with £ and whether 
or not no spreading and catastrophic self-focusing of the anomalous pulse develop. Nevertheless, 
the currently available numerical results suggest that a self-trapped solution can exist in the 
configuration under discussion. Further calculations should clarify this question. 

Note that a configuration of two simultaneously propagating pulses could also be used in 
optical compression techniques since, as is seen from figure ||(a), for some particular values of the 
longitudinal distance £ the temporal width of the anomalous pulse decreases by about five times 
the initial width. 

6 Conclusions 

In this paper properties of two pulses propagating simultaneously in different dispersion regimes, 
i.e. anomalous and normal, in a Kerr-type planar waveguide are considered. The propagation 
is described by two coupled NSEs. The interaction between pulses is assumed to be limited to 
cross-phase modulation. Four wave mixing is neglected, i.e. no energy transfer between pulses is 
taken into account. The accuracy of another assumption used in the analysis, the omitting of the 
difference of group velocities of the pulses, is discussed in appendix B. Our analysis is based on 
the variational method and numerical simulations. 

First we studied the influence of the parameters of the pulse propagating in a normal dispersion 
regime on the threshold of catastrophic self-focusing of the pulse with an anomalous dispersion. 
We observed that catastrophic self-focusing of the pulse propagating in an anomalous dispersion 
regime can be arrested by the pulse propagating in a normal dispersion regime when the strength 
of nonlinearity is sufficiently large, K2 > nxiow and the dispersion-to-diffraction ratio satisfies the 
relation: \<Jxiow{n2)\ < {02] < Wxupp{n2)\- In this notation X = V(X = N) concerns the results 
obtained in the variational method (numerical simulations). We also investigated whether the 
nonlinear coupling between pulses can cause catastrophic sclf-focusing of the pulse propagating in 
a normal dispersion regime. The variational method indicates that when catastrophic self-focusing 
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Figure 5: The results obtained using the variational method displaying the dependence of the 
temporal, w T , and spatial, w^, widths of the anomalous pulse (a) co-propagating with the normal 
pulse (b), for the following parammeters of the system: K\ = 1.0, o\ = 1.0, «2 = 2.0, a 2 = 0.0. 
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Figure 6: The results obtained using the numerical simulations displaying the dependence of the 
temporal, w T , and spatial, w^, widths of the anomalous pulse (a) co-propagating with the normal 
pulse (b) for the following parammeters of the system: m = 1.0, u\ = 1.0, k 2 = 2.0, a 2 = 0.0. 
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of the anomalous pulse occurs, the normal pulse can display, depending on the parameters of the 
system, two different characteristics: (i) both widths of the pulse initially decrease reaching a 
minimum on a certain distance of propagation and then they start to increase, (ii) the spatial 
width of the pulse vanishes to zero on a finite distance of propagation whereas the temporal width 
initially decreases reaching a minimum on a certain distance of propagation and then it increases. 
The occurence of catastrophic self-focusing of the normal pulse has not been observed in the 
numerical simulations. Therefore, we can conclude, based on the variational method and the 
numerical simulations, that catastrophic self-focusing of the anomalous pulse does not necessarily 
lead to catastrophic self-focusing of the normal pulse. 

We found also, using the numerical simulations, that the presence of the pulse propagating 
in a normal dispersion regime can lead to spatio-temporal splitting of the pulse propagating in 
an anomalous dispersion regime. Recall that splitting of an anomalous pulse into several pulses 
does not occur when it propagates as a single pulse. 

Finally, we considered the limiting case of vanishing dispersion of the pulse propagating in a 
normal dispersion regime with parameters of the pulses chosen in such a way that catastrophic 
self-focusing does not occur, i.e. that the conditions «i < Kycat = 1 + ffi (in the variational 
method) and K\ < KNcat ~ 0.885 + o\ (in the numerical simulations) are satisfied. The main 
motivation was to see whether such a configuration can lead to a stable self-trapped propagation 
of a pulse with anomalous dispersion. The positive answer was obtained within the variational 
method which confirms that neither spatio-temporal spreading nor catastrophic self-focusing of 
the anomalous pulse can develop thus giving rise to a self-trapped solution. Note that this kind 
of stabilization is similar to that which has been found earlier in media with saturation-type 
nonlinearity Although the existing data supports the existence of a self-trapped solution, 

conclusive results require labourious simulations at high values of the longitudinal variable £ and 
are not yet available (work in progress). 

Note, in conclusion, that the existence of a stable self-trapped solution could be useful, for 
example, in optical switching devices. The configuration of two simultaneously propagating pulses 
in a planar waveguide could also be of use in optical compression techniques. 
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Appendix A 

The notation in equations (|l|) and (||a,b) is as follows 0: ( = zjzuFi is the longitudinal co- 
ordinate normalized to the Fresnel diffraction length of the anomalous pulse, £ = x/w\ is the 
spatial transverse coordinate normalized to the initial spatial width of the anomalous pulse, 
t = (t — 0i z)jt\ is the local time normalized to the initial temporal width of the anomalous 
pulse. The parameters (Xj — Zdfi/ ZDSj, A* = zdfi/ 'zdf2, t = A1/A2 = 1^2/^1 denote, re- 
spectively, the dispersion-to-diffraction ratio, the ratio of the Fresnel diffraction length of the 
anomalous pulse to the Fresnel diffraction length of the normal pulse and, finally, the ratio 
of the carrier frequency of the anomalous pulse to the carrier frequency of the normal pulse. 

: = y/KjUj(£,T,£)/Ujo denotes the normalized amplitude of the j-s pulse, where Uj(£,T,£) 
is the amplitude of the slowly varying envelope of the electric field, Ujo := 1^ (0,0,0) is the di- 
mensionless initial peak amplitude. The parameter Kj := {zdfj / znlj) 2 defined as the strength 
of nonlinearity of j-s pulse is proportional to the nonlinear part, 722, of the refractive index of 
a medium, n :— hq + n2\Uj\ 2 , to the initial peak intensity, | C/^ 1 2 7 and to the square of the 
spatial width of the pulse, w%. Note also that in the case of the (l-fl)-dimensional NSE, 
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i.e. when Oj = 0, the quantity y/77J can be interpreted as the order of a spatial soliton, so 
that a first-order soliton arises when kj — 1 |6l| . The dispersive terms are defined as fol- 
lows: f3j QS> := f3^(ujj) — LOj/c is the wavenumber, f3j := d(3jdui\ U j =u}j = l/v g j is the reverse 
group velocity, and (3^ := d 2 (3 / doJ 2 \ UJ=U j j is the group velocity dispersion. The parameters 

zdfj ■= /3] 0) n (^)wf , z D Sj '■= \ z NLj ■= w 3V / nQ/(2n 2 \U 0j \ 2 ), Wj, tj denote, respectively, 

the Fresnel diffraction length, the dispersive length, the nonlinear length, the initial spatial width 
and the initial temporal width of the j-s pulse. In the above notation j = 1,2, where the subscript 
j = 1 (j = 2) refers to the anomalous (normal) pulse. 

Appendix B 

Since we have assumed that pulses have different wavelengths and different group velocity dis- 
persions, it is physically evident that they should also have different group velocities. Therefore, 
the assumption that the difference of the group velocities of the pulses vanishes is a simplification 
accepted in this paper and should be viewed as a first step of the analysis. When this difference 
does not vanish the pulses propagate with different velocities and the overlap between them 
decreases with the longitudinal variable. Therefore, the nonlinear coupling between them also 
decreases. In the limiting case of the difference of the group velocities of the pulses approaching 
infinity the coupling between pulses becomes zero and the problem of simultaneous propagation 
of two pulses reduces to the case when they propagate separately. 

However, we believe that the inclusion of a small difference of group velocities of the pulses, 
which should be studied numerically, will not cause qualitative changes in the results of this paper, 
such as the possibility of an arresting of catastrophic self-focusing of the pulse propagating in 
an anomalous dispersion regime by the influence of the pulse propagating in a normal dispersion 
regime. The only difference we expect is a change of the values of the parameters, ONiowiVNupp, 
i^Nuppi which describe the threshold of catastrophic self-focusing for fixed values of a\ and K\. 
These quantitative changes would be proportional to the value of the difference of the group 
velocities of the pulses. 

References 

[1] Y. Silberberg, Opt. Lett. 15, 1282 (1990). 

[2] A. T. Ryan and G. P. Agrawal, Opt. Lett. 20, 306 (1995). 

[3] M. Desaix, D. Anderson, and M. Lisak, J.Opt.Soc.Am. B. 8, 2082 (1991). 

[4] V. E. Zakharov and A. B. Shabat, Zh. Exp. Teor. Fiz. 61, 118 (1971). 

[5] G. Fibich, Phys. Rev. Lett. 76, 4356 (1996). 

[6] S. N. Vlasov, V. A. Petrishchev, and V. I. Talanov, Izv. Vuz Radiofiz. 14, 1353 (1971). 

[7] J. J. Rasmussen and K. Rypdal, Physica Scrypta 33, 481 (1986). 

[8] X. D. Cao, G. P. Agrawal, and C. J. McKinstrie, Phys. Rev. A 49, 4085 (1994). 

[9] J. T. Manassah, P. L. Baldeck, and R. Alfano, Opt. Lett. 13, 1090 (1988). 

[10] P. Chernev and V. Petrov, Opt. Comm. 87, 28 (1992). 

[11] M. Pietrzyk, Opt. and Q.E. 29, 579 (1997). 

[12] M. Karlsson, Phys. Rev. A 46, 2726 (1992). 

[13] D. Suter and T. Blasbberg, Phhys. Rev A 47, 256 (1993). 

[14] A. L. Dyshko, V. N. Lugovi, and A. M. Prokhorov, Sov. Phys. JETP 34, 1235 (1972). 

[15] E. Yablonovitch and B. Bloembergen, Phys. Rev. Lett. 29, 907 (1972). 

[16] S. Henz and J. Hermann, Phys. Rev. E 53, 4092 (1996). 

[17] V. I. Karpman, Phys. Rev. E 53, R1336 (1996). 

[18] S. Blair and K. Wagner, Opt. and Q.E. 30, 697 (1998). 

[19] N. Akhmediev, A. Ankiewicz, and J. M. Soto-Crespo, Opt. Lett. 18, 411 (1993). 

[20] G. Fibich and G. C. Papanicolaou, Opt. Lett. 22, 1379 (1997). 

[21] I. Gurwich, J. Phys. D 30, 2183 (1997). 



17 



[22] M. D. Feit and J. J. A. Fleck, J. Opt. Soc. Am. B 5, 633 (1988). 

[23] A. P. Sheppard and M. Healtermand, Opt. Lett. 23, 1820 (1998). 

[24] F. H. Berkshire and J. D. Gibbon, Stud. Applied Math. 69, 229 (1983). 

[25] A. G. Litvak, T. A. Petrova, A. M. Sergeev, and A. D. Yunakovskii, Sov. J. Plasma Phys. 9, 287 
(1993). 

[26] M. Pietrzyk, PhD thesis, Warsaw 1998 (in Polish, unpublished) . 
[27] B. Gross and J. T. Manassah, Opt. Comm. 126, 269 (1996). 

[28] D. Burak and R. Binder, in Proc. Quantum Electronics and Laser Science Conference (OSA, USA, 
1996). 

[29] J. E. Rothenberg, Opt. Lett. 17, 1340 (1992). 

[30] P. Chernev and V. Petrov, Opt. Lett. 17, 172 (1992). 

[31] G. G. Luther, J. V. Moloney, A. C. Newell, and E. M. Wright, Opt. Lett. 19, 862 (1994). 

[32] J. E. Rothenberg, Opt. Lett. 17, 583 (1992). 

[33] M. Trippenbach and Y. B. Band, Phys. Rev. A 56, 4242 (1997). 

[34] A. A. Zozulya, S. A. Diddams, A. G. VanEugen, and T. S. Clement, Phys. Rev. Lett. 82, 1430 
(1999). 

[35] A. A. Zozulya and S. A. Diddams, Opt. Express 4, 336 (1999). 
[36] J. K. Ranka and A. L. Gaeta, Opt. Lett. 23, 534 (1998). 

[37] J. K. Ranka, R. W. Schirmer, and A. L. Gaeta, Phys. Rev. Lett. 77, 3783 (1996). 

[38] S. A. Diddams, H. K. Eaton, A. A. Zozulya, and T. S. Clement, Opt. Lett. 23, 379 (1998). 

[39] C. J. McKinstrie and D. A. Russell, Phys. Rev. Lett. 61, 2929 (1988). 

[40] L. Berge, Phys. Rev. E 58, 6606 (1998). 

[41] A. S. Desyatnikov and A. I. Maimistov, J. Exp. Theor. Phys. 86, 1101 (1998). 

[42] O. Bang, L. Berge, and J. J. Rasmussen, Phys. Rev. E 59, 4600 (1999). 

[43] G. P. Agrawal, Phys. Rev. Lett. 64, 2487 (1990). 

[44] K. Konno and H. Suzuki, Phys. Scripta 20, 382 (1979). 

[45] R. McLeod, K. Wagner, and S. Blair, Phys. Rev. A 52, 3254 (1995). 

[46] F. Vidal and T. W. Johnston, Phys. Rev. Lett. 77, 1282 (1996). 

[47] S. Gatz and J. Herrmann, J. Opt. Soc. Am. B 14, 1795 (1997). 

[48] V. Skarka, V. I. Berezhiani, and R. Miklaszewski, Phys. Rev. E 56, 1080 (1997). 

[49] M. Segev, B. Crosignani, and A. Yariv, Phys. Rev. Lett. 68, 923 (1992). 

[50] M. F. Shih and M. Segev, Phys. Rev. Lett. 78, 2551 (1997). 

[51] B. A. Malomed et al, Phys. Rev. E 56, 4725 (1997). 

[52] L. Berge, V. K. Mezentsev, J. J. Rasmussen, and J. Wyller, Phys. Rev. A 52, R28 (1995). 
[53] K. Hayata and M. Koshiba, Phys. Rev. Lett. 71, 3275 (1993). 

[54] D. Mihalache, D. Mazilu, L. C. Crasovan, and L. Torner, Opt. Comm. 137, 113 (1997). 

[55] P. Agin and G. I. Stegeman, JOSA B 14, 3162 (1997). 

[56] W. E. Torruellas et al, Phys. Rev. Lett. 74, 5036 (1995). 

[57] X. Liu, L. J. Qian, and F. W. Wise, Phys. Rev. Lett. 82, 4631 (1999). 

[58] L. Berge, O. Bang, J. J. Rasmussen, and V. K. Mezentsev, Phys. Rev. E 55, 3555 (1997). 

[59] S. Darmanyan, I. Relke, and F. Lederer, Phys. Rev. E 55, 7662 (1997). 

[60] D. Anderson and M. Bonnedal, Phys. Fluids 22, 105 (1979). 

[61] G. P. Agrawal, "Nonlinear Fiber Optics" (Academic Press, London, 1989). 

[62] D. Anderson, M. Lisak, and T. Reicher, J. Opt. Soc. Am. B 5, 207 (1988). 

[63] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, "Numerical Recipes in Fortran" 

(Cambridge University Press, Cambridge, 1992). 
[64] F. Cornolti, M. Lucchesi, and B. Zambon, Opt. Comm 75, 129 (1990). 
[65] M. Desaix, D. Anderson, and M. Lisak, Phys. Rev. E 50, 2253 (1994). 
[66] L. Berge et al, J. Opt. Soc. Am. B 13, 1879 (1996). 
[67] D. Burak and W. Nasalski, J. Tech. Phys. 36, 199 (1995). 
[68] D. Anderson, Phys. Rev. A 27, 3135 (1983). 
[69] D. Burak, Opt. Appl. 21, 3 (1991). 



18 



